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The critical behavior of quantum Hall transitions in two-dimensional disordered electronic systems 
can be described by a class of complicated, non-unit ary conformal field theories with logarithmic 
correlations. The nature and the physical origin of these logarithmic correlation functions remain 
however mysterious. Using the replica trick and the underlying symmetries of these quantum critical 
points, we show here how to construct non-perturbatively disorder-averaged observables in terms 
of Green’s functions that scale logarithmically at criticality. In the case of the spin quantum Hall 
transition, which may occur in disordered superconductors with spin-rotation symmetry and broken 
time reversal invariance, we argue that our results are compatible with an alternative approach 
based on supersymmetry. The generalization to the Integer quantum Hall plateau transition is also 
discussed. 


I. INTRODUCTION 

Random impurities in systems of non-interacting elec¬ 
trons can induce a transition between metallic (delocal¬ 
ized) and insulating (localized) phased. More than fifty 
years after its discovery, the field of Anderson localization 
remains very activ^. Even though the standard scaling 
theory of localizatioiP predicts the absence of extended 
states in disordered non-interacting electronic systems in 
two dimensions, a well-known exception to this rule is 
provided by the transition between plateaux in the In¬ 
teger Quantum Hall Effect (IQHE), a quantum critical 
phenomena that was predicted theoretically and observed 
experimentally a few decades agcP^^. 

The physics of non-interacting electrons moving on 
a two-dimensional plane with a perpendicular magnetic 
field and a random potential can be described by a 
two-dimensional sigma model on the manifold = 
U{2n)/U{n) x U{n) in the limit n ^ (Pl. The sigma 
model can be argued to be asymptotically free, and it 
flows to strong coupling, suggesting complete localiza¬ 
tion in dimension d = 2 (this is related to the absence of 
Goldstone phase in d = 2). The reason why the IQHE 
transition is allowed is well understood in terms of an 
additional topological 0 term - that can be written be¬ 
cause 7r2(M^) = Z - in the sigma model that induces 
a transition at ^ = tt. It is by now well-admitted that 
the quantum critical point corresponding to this IQHE 
transition - the strong coupling fixed point of the sigma 
model nt 0 = IT - should be effectively described by a 
Conformal Eield Theory (GET) in l+I dimensions (see 
e.g. Refs. I8I9|) . 

This effective description of a disordered electronic sys¬ 
tem in 2-|-ld in terms of a GET in I-fld comes at the 
price of losing unitarity, and as a result, the critical 
properties of the IQHE transition remain quite poorly 
understood: the critical exponents are known only ei¬ 
ther from experiments or numerical simulations, and a 
low energy field theory description of the critical point is 
still missing. Because of the underlying replica or super- 
symmetric description of the disordered critical point, it 


is however easy to ar gue t hat the central charge of this 
GET should be c = Cpniiil _ because the partition func¬ 
tion of the problem is essentially trivial Z = 1 - and as 
such should be a Logarithmic Gonformal Eield Theory 
(LGET)P^lli^ (the only unitary non-logarithmic GET at 
c = 0 being trivial). LGETs are characterized by the 
non-diagonalizability of the scale transformation genera¬ 
tor, corresponding to indecomposable (reducible but not 
fully reducible), non-unitary representations of the Vi- 
rasoro algebra. This property leads to the appearance 
of logarithms in correlation functions and to “indecom¬ 
posable” operator product expansions. These quantum 
field theories are daunting for many reasons: they are 
necessarily non-unitary, they are typically non-rational 
{i.e. they involve infinitely many primary fields), and 
from an algebraic point of view they require going be¬ 
yond the description in terms of irreducible representa¬ 
tions familiar to physicists, and are instead built out of 
large, complicated indecomposable representations that 
are very hard to classify mathematicall^ff^. Most tools fa¬ 
miliar to physicists in this context also have to be recon¬ 
sidered: for instance the Mermin-Wagner theorem does 
not apply, and conformal invariance and global group 
symmetry does not imply Kac Moody symmetry. Eor 
example, whereas the 0(3) sigma modeP^ on the coset 
S‘^ = SU{2)/U{1) at topological angle 0 = TT flows to a 
Wess-Zumino-Witten (WZW) GET with enlarged sym¬ 
metry SU{2)i that is “easily” manageabl^^^, the GET 
describing the IQHE transition (also obtained from a 
seemingly similar sigma model) is most likely not a WZW 
model. Despite the recent progress in the understanding 
of such LGETs both from lattice modelJ^^M^ and more 
abstract algebraic approaches (see e.g. Refs. [2T]-[26|), all 
examples of LGETs that are slowly getting under control 
are still pretty far from the expected complexity of the 
LGET describing the IQHE plateau transition. 

Even if a detailed understanding of this logarithmic 
GET seems unfortunately out of reach for now, it should 
be possible to understand what physical observables make 
the theory logarithmic. One of the main features of a 
LGET is indeed the existence of logarithmic correlations 
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at criticality - where one usually expects only power- 
laws; it is therefore very natural to ask what kind of 
disorder-averaged observables in the IQHE show loga¬ 
rithmic correlations at the plateau tra nsitio n. A gen¬ 
eral mechanism was proposed by Card^W^^ to explain 
the emergence of logarithmic correlations in disordered 
systems from the replica trick, but this argument has 
to be adapted to the more complicated symmetry struc¬ 
ture of the IQHE critical point. This approach was re¬ 
markably successful to understand the physical origin of 
logarithmic correlati ons in “simpler” problems including 
the 0{n 0) modeP^^, or percolatioiP^ for example. 

Even though the IQHE transition is admittedly much 
more involved, it is still natural to expect the (heuris¬ 
tic) arguments of Refs. I13I27I to shed some light on the 
physical origins of logarithmic correlations in Anderson 
transitions. 

Eollowing Cardy, our general strategy in this paper will 
be to use the replica trick to express disorder-averaged 
observables (Green’s functions) in terms of correlators of 
n species of fermionic operators, with n the number of 
replicas, in a theory with a global group symmetry G{n) 
- typically the unitary group U{n). Using the represen¬ 
tation theory of this global symmetry group, we classify 
the different lattice observables and deduce the general 
form of their correlation functions using symmetry ar¬ 
guments only. The critical properties of the physical, 
disordered system can then be understood, in principle, 
from the formal n ^ 0 limit of the critical point (if it 
exists) of these replicated theories for generic n integer. 
The resulting correlation functions usually have 1/n poles 
and are therefore ill-defined in the physical limit n ^ 0. 
These apparent divergences in the limit n ^ 0 can be 
cured but they lead to logarithms in the limit. The ba¬ 
sic mechanism underlying the emergence of logarithms in 
the replica limit is quite easily understood: starting from 
two power-law functions /n and with 

coinciding critical exponents for n = 0, both diverging as 
1/n in the limit n ^ 0, a well-defined quantity can be 
obtained as 


^-2Ai(n) _ ^-2A2(n) 

lim - 

n^O n 


= hvr 


( 1 ) 


where n = 2 


is a universal number. We 


n=0 


will say that the two scaling operators corresponding to 
the critical exponents Ai(n) and A 2 (n) are “mixed” at 
n = 0, the precise sense of this statement is that the scale 
transformation operator becomes non-diagonalizable at 
n = (P. These correlation functions with logarithmic 
terms can in turn be argued to correspond to logarith¬ 
mic operators in a LCET, so that scale invariance is fully 
preserved at the critical point despite the presence of non- 
algebraic correlations. The physical meaning of these log¬ 
arithmic correlations can then be inferred using Wick’s 
theorem, thus yielding disorder-averaged observables in 
terms of Green’s functions that should behave logarith¬ 
mically at criticality. 


The apparent simplicity of this program is somewhat 
deceiving: the crucial point is to identify properly the 
symmetry G{n) of the critical point for generic integer n. 
Note in particular that our approach is drastically differ¬ 
ent from the classification of operators based on the sigma 
model formulation of the transitiorP^^^ (see also|9]): af¬ 
ter all, the theory describing the critical point probably 
has nothing to do with Pruisken’s sigma modeP (or its 
supersymmetric variantJ^ - cf. the above discussion of 
the WZW model that arises in the strong limit coupling 
of the 0(3) sigma model. On the other hand, the actual 
symmetry of the critical point is probably much larger 
that the naive G{n) symmetry identified from the action 
of the replicated theory. This means that two operators 
that transform under different irreducible representations 
of G{n) (and that could eventually be “mixed” in the 
limit n ^ 0) could actually be part of the same multiplet 
of fields for the actual, larger symmetry. Erom a prac¬ 
tical point of view, this means that the two exponents 
Ai and A 2 in eq. 0 could be identical, implying that 
/i: = 0, that is, no logarithmic term. Despite this impor¬ 
tant issue, similar approaches applied to simpler cr itical 
points such as the 0{n 0) model or percolatiorP^^^ 

(obtained as the Q ^ 1 limit of the Q-state Potts model 
with perturbation group Sq symmetry) were remarkably 
successful to classify the operator content and the loga¬ 
rithmic correlations of these critical points, even though 
it is known that the actual symmetr}!^ of these theo¬ 
ries is in fact much larger than 0{n) or Sq. Even though 
quantum Hall transitions are admittedly much more com¬ 
plicated than say the 0{n = 0) model or percolation, 
it is natural to try to extend this approach to critical 
points that are also much more exciting from a physi¬ 
cal perspective. In the following, we will see that such 
a symmetry-based analysis with the simplest symmetry 
groups possible yields very sensible results, in agreement 
with the existing literature on the subject. In the case of 
the Spin Quantum Hall Effect (SQHE), a close cousin of 
the Integer Quantum Hall transition where SU (2) spins 
replace U (1) charges, we predict the existence of logarith¬ 
mic corrections in some observable that is found to be in 
agreement with another approach, based on supersym¬ 
metry and a mapping onto percolatiorP^. Eor the IQHE 
transition, we show that the disorder average of a sim¬ 
ple combination of Green’s functions should scale exactly 
as ^ logr at the critical point, without any power-law 
contribution; a concrete prediction that could be veri¬ 
fied numerically using the Ghalker-Goddington network 
modeP^. 

We emphasize that extrapolating the limit n ^ 0 of 
the replica limit is notoriously complicated in the context 
of the IQHE - a fact that can essentially be traced back 
to the fact that the transition at ^ = tt is believed to be 
of first order for n > 2, making very h ard to extrapolate 
reliably physical results to n = CPHEU Nevertheless, the 
replica trick seems more reliable when it comes to predict¬ 
ing the existence of logarithmic correlations (compared to 
say, computing the exact value of universal quantities like 















3 


critical exponents in the limit n ^ 0). The reason for this 
is that the basic mechanism that predicts the emergence 
of logarithms for n ^ 0, eq. 0 - does not rely on the 
precise value of the critical exponents or on details of the 
extrapolation, so that one can essentially ignore that the 
transition becomes of first order for larger n. This is well 
illustrated in the simple case of the n ^ 1 limit of the 
sigma model with topological ^-term, which also 
has a phase transition at 6> = tt. Despite the fact that the 
transition is also believed to be of first order for n > 2, 
we will see below that the replica trick based on the (at 
least) U (n) symmetry of the critical point leads to results 
that are in agreement with a more rigorous approaches 
that are available in this cas^^^^^^EH. 

The remainder of this paper is organized as follows. 
In section [n| we analyze the n ^ 1 limit of the 
sigma model, a quantum field theory describing the crit¬ 
ical behavior of dense loops with fugacity n, a problem 
apparently unrelated to Anderson transitions. We de¬ 
scribe in details how logarithmic correlations arise in 
the n ^ 1 limit (corresponding to the classical perco¬ 
lation problem). Building on these results, we analyze 
the replica limit of the SQHE transition in section |III| 
and predict the existence of logarithmic correlations at 
the critical point. We show that the same result can also 
be obtained from the analysis of the percolation problem 
in Sec. [H] (see also ref. EH using the supersymmetry ap¬ 
proach developed in Ref. [36| We then extend the replica 
approach to the IQHE case in section Eg and construct 
an observable that should behave purely logarithmically 
at the critical point. Einally, section |V| contains a dis¬ 
cussion of the results and concluding remarks. The main 
results of our analysis are given by eqs. (29) and (31) for 
the SQHE, and eqs. (38) and (39) for the IQHE - see 
also eq. 


II. WARM-UP: n 1 LIMIT OF THE CP^”^ 
SIGMA MODEL, PERCOLATION AND DENSE 
LOOP MODELS 


We begin with a discussion of logarithmic correlations 
in a model of dense loops with SU (n) symmetry, a prob¬ 
lem apparently completely unrelated to Anderson tran¬ 
sitions. However, most results of this paper regarding 
Anderson transitions rely heavily on this section. The 
aim of this first section is three-fold: (1) it will provide 
us with a simple example where the emergence of loga¬ 
rithmic correlations can be understood in details, (2) it 
introduces the basis of SU (n) representation theory that 
will be crucial in the following, and (3) it proves the exis¬ 
tence of a logarithmic observable in a critical dense loop 
model with fugacity n = 1 (percolation), a model that 
was showiP^ to describe a certain class of observables in 
the SQHE transition using the so-called supersymmetry 
trick. We will show in the next section how this result 
can be interpreted in the SQHE language using this su¬ 
persymmetry mapping. 


A. Dense loop models and SU(n) symmetry 

Let us consider a dense loop model on the square lat¬ 
tice, with fugacity (Boltzmann weight per closed loop) 
n G M. An example of configuration of such a loop model 
is shown in Eig. where the dense loops are drawn in 
blue. These loops can be thought of as “hulls” of per¬ 
colation clusters shown in red in Eig. more precisely, 
the red clusters are percolation clusters only for n = 1; 
for generic n, they coincide with the so-called Eortuin- 
Kasteleyn clusters that arise in the high-temperature ex¬ 
pansion of the Potts modeP^. There are two types of 
plaquettes in Eig. we will denote the plaquettes corre¬ 
sponding to the propagation of two blue lines along the 
vertical (imaginary time) direction (with a vertical red 
line) by 1 , and we will soon see why we use this notation 
of identity operator. On the other hand, the other type 
of plaquette with a horizontal red line are denoted by an 
operator e that contracts two blue lines and creates a new 
pair of lines in imaginary time. We will see in the follow¬ 
ing that it is useful to think of these loops as worldlines 
of bosonic particles that get created and annihilated by 
the operators e during the imaginary time evolution. 

Labeling by i the columns in Eig.[^ the partition func¬ 
tion is computed as follows: if i is even, plaquettes of type 
1 are weighted by (1 — p^) whereas plaquettes associated 
with the operator e get a weight If ^ is odd, the 

same applies with the replacement pa 1 — Pb- Einally, 
closed loops carry a Boltzmann weight n. The isotropic 
line corresponds to p^ = Pb and one can show that the 
system is critical for p^ = 1 — Pb and —2 < n < 2 (for 
other values of n, there is a first order transition). More 
formally, we introduce the transfer matrix 

T= Yl{{l-pA)'i-i+PAei) P ipB'i-i + {l-PB)ei), (2) 

i odd i even 

which constructs two rows of the loop model. A full 
loop configuration such as the one in Eig. can then 
be constructed by successive iterations of this transfer 
matrix. In the strongly anisotropic limit Pa ^ 0 with 
Pa/(1 — Pb) fixed, the transfer matrix can be recast as 
T exp(—y^pA(l —Pb)H)^ with the (H-l)D quantum 
Hamiltonian 


H = -e^ Ci- e ^ y] Ci, (3) 

i odd i even 


with e = \JpaI (1 — Pb) = 1 at the critical point. 

We now discuss a SU{n) representation of this dense 
loo p mode l for integer n. There are several ways to do 
thiJlSIlSIlll follow here Ref. |40l On each site i (to 
simplify the notations, we implicitly work on a given 
imaginary time slice of the system so that i labels the 
columns of the network as in eqs 0 and ([^), we in¬ 
troduce n bosonic operators with commutation 

relations [b^\b^j] = Sij6^^ and a,b = l,...,n. These 
operators satisfy the following constraint = 1 
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FIG. 1: Configuration of a dense loop model with fugacity 
n = 1, and corresponding percolation clusters in red. 

(one particle per site), so that the n states \a) = |0), 

with |0) the Fock space vacuum, form the fundamental 
(resp. anti fundamental or dual) representation of the 
Lie algebra su(n) for i even (resp. odd), with the gen¬ 
erators acting as Qf^ = (resp. Qf^ = 

The operator is then defined as the projector onto 
the singlet in the tensor product between fundamental 
and dual if i is even (dual and fundamental for i odd) 
Ci = — that has the familiar form of a 

Heisenberg-like coupling. It is then straightforward to 
verify that the partition function computed as the trace 
of powers of the transfer matrix (© admits a graphical 
expansion in terms of non-intersecting dense loops - the 
worldlines of the bosons, with the weight of each closed 
loop given by the trace of the identity operator in the 
fundamental representation (number of particles flowing 
around the loop), that is, n. Therefore, the transfer ma¬ 
trix (or the corresponding quantum Hamiltonian) of our 
loop model can naturally be endowed with a SU (n) sym¬ 
metry (when n > 1 is integer). 

The quantum field theory description of these loop 
models (or of the corresponding SU{n) quantum spin 
chains) is then provided by a sigma modeP^^^^^ on the 
coset = U{n)/{U{1) x U{n — 1)) with topological 

angle 0. We introduce n bosonic fields Za^ a = 1,... ,n 
subject to the constraint z^Za = 1 (implicit summa¬ 
tion over repeated indices implied), modulo U{1) phases 
SO that z G With 

this field content, the Lagrangian density of the sigma 
model reads 

1 iO 

^ ~ 2^ T ^ (4) 

where is a U{1) gauge field and = 9^ + ia^. 
Integrating out the abelian gauge field yields = 
I — {d^z\)za). The ^-angle contribution is a topo¬ 

logical term associated with the non-trivial homotopy 


group 7r2(CP’^“^) = Z. 

If we now think of n as a real parameter after a naive 
analytical continuation, it can be shown from the p func¬ 
tion of the sigma modeP^that the coupling flows to larger 
values at larger length scales (lower energies) for n > 0. If 
^ 7 ^ 0, the system flows to a massive phase with restored 
U(n) symmetry, but there is a transition at topological 
angle 0 = tt. For n < 2, this transition is believed to be of 
second order and should therefore be described by a con¬ 
formal field theory (CFT). Many properties of this CFT 
as a function of n are knowiP^, including for example the 
expression of the central charge and of critical exponents, 
and in the following we will argue that logarithmic corre¬ 
lations arise naturally from the analytic continuation of 
the SU (n) symmetry at the value n = 1 (corresponding 
to the percolation problem). Note that although we will 
focus in what follows on this second quantized descrip¬ 
tion of SU (n) loop models in two dimensions, most of our 
results also apply to higher dimensional loop modelP^^^. 

B. Observables in SU{n) loop models 

In order to classify the observables of the loop mod¬ 
els, we use the underlying SU{n) symmetry of the crit¬ 
ical point. Note that the actual symmetry at the criti¬ 
cal point might be much larger than SU(n) - for these 
simple loop models, the symmetry is knowiP^ and is in¬ 
deed much larger than SU{n) . As a result, observables 
transforming under different irreducible representations 
of SU (n) might be part of a larger multiplet - and in par¬ 
ticular would have the same scaling dimension - under 
the eventual larger symmetry. We will nevertheless keep 
working with the SU(n) symmetry and analyze carefully 
the results in the end. We will focus for concreteness on 
the two-dimensional second quantized model introduced 
above, but almost id entica l results apply to higher di¬ 
mensional loop modelJ^^Ell. 


J. Observables acting on N — 1 site 

Let us start by classifying the hermitian opera¬ 
tors acting on a single edge of the form 0{i) = 

6 ^a6^1(^)^6(^)- Clearly, the total number of bosons 
is invariant under SU{n), and form 
an irreducible unidimensional representation. Within the 
context of our loop models, 4>(i) = 1 is the identity op¬ 
erator, with scaling dimension A<^ = 0. The remaining 
- 1 operators (f>ab{i) = bl{i)bb{i) - ^ 
transform irreducibly and form the adjoint representation 
of SU (n). They are traceless in the sense that they satisfy 
'^ab ^ab^^abi'^) = 0- This decomposition can be written in 
terms of Young diagrams as [1] (8) [1]"*" = [1] (8) [n — 1] = 
[n] 0 [n—1,1]. Here and in the sequel. Young diagrams are 
denoted as [Ai, A 2 ,... ], where A^ is the number of boxes 
in the column. Here, [1] is the fundamental (defin¬ 
ing) representation of dimension n and [1]* its conjugate 
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(dual), and [n — 1,1] is the adjoint. Using the fact that 
tT(j)ah{i) = 0 , where the trace is taken in the fundamental 
representation of SU{n) generated by \a) = bl^\0)^ it is 
then straightforward to show that the two-point function 
of (j)ab can be expressed as 

{4>aa{i)4>bb{j)) = - (Sab “ “) ( 5 ) 

n \ n J 

where is the probability than i and j belong to 

the same loop. The factor (^6ab ~ is completely fixed 
by the condition (paa = 0- Note that we have fixed 
some indices for simplicity but the full correlation func¬ 
tion ( 0 a 6 (^) 0 cd(i)) can of course be computed similarly. 
One can also readily check that {^{i)(l)ab{j)) = 0 , and 
= (^) = 0 . 

This exact lattice formula makes the nature of the op¬ 
erator 0 very transparent: it creates a propagating loop, 
or more precisely, it creates two “legs” at site i, one in¬ 
coming and one outgoing. In two dimensions, this observ¬ 
able is known as a “ 2 -leg watermelon” operator for this 
reason. Let A 0 (n) be its scaling dimension, known ex¬ 
actly in two dimensionP^ A 0 (n) = 1 — with g ^ [ 2 ,4] 
given by = 2(1 + cos ^). In the continuum limit, we 
expect the following scaling 

{<i>aa{n)4>bb{rj)) = A{n) ( 6 ) 

where A{n) is some non-zero function of n, and r = |r^ — 

rj\- 

To summarize, the SU{n) symmetry allowed us to 
identify two distinct scaling operators, the identity and 
a 2-leg watermelon operator. While the representation 
theory analysis was performed for n integer, the correla¬ 
tion functions (HI,® make sense for generic n. Although 
nothing particularly exciting happens in the limit n ^ 1 , 
one can already see that the limit n ^ 0 will be singular 
from eq. Ill- As we will see more explicitly in the fol¬ 
lowing, poles like the one in eq. in lead to logarithmic 
correlations. In this case, one can check that the scaling 
dimension A 0 (n) vanishes as n ^ 0 , indicating a “mix¬ 
ing” of the operator (j)ab with 4> (the identity) at n = 0 . 
This mixing can also be traced back to the dimension 
— 1 of the representation of (j) that becomes formally 
— 1 in the limit n ^ 0. We will not describe this limit 
n ^ 0 here, but instead turn to more complicated opera¬ 
tors that will have singular limits as n ^ 1. We will then 
show more explicitly how these poles lead to logarithmic 
correlations. 


2. Observables acting on N — 2 sites 

We now turn to observables acting on two nearest 
neighbor edges i and i'. For simplicity, we consider a 
coarse-grained picture in which i and i' are two sites near 
position Vi (see Fig. |^, corresponding to the fundamen¬ 
tal representation of SU{n) (even columns). We wish 
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FIG. 2: Schematic representation of the coarse grained two- 
point function (0{ri)0{rj)) of observables acting on A = 2 
sites. The two sites i and i' (resp. j and j') are chosen in the 
infinitesimal neighborhood of the position n {resp. rj). 

to consider operators of the type b\{i)bi){i)bl{i')bd{i'). 
It is natural to enforce a symmetry between i and i', 
and we further restrict our study to operators that are 
symmetric under the exchange (a, 6) (c, d), and that 

vanish if a = c or 6 = d. The last two conditions are 
conveniently satisfied by operators that are antisymmet¬ 
ric under the exchanges a ^ c or b ^ d. We there- 
fore define Tabcd{ri) = -Qab{i)Qcd{i') + Qad{i)Qcb{i') - 
Qab{i')Qcd{i) + Qad{i')Qcb{i), with Qab{i) = bl{i)bb{i). 
These operators transform under SU{n) as [2] (8) [2]* C 
([1] (g) [1]*)*^^, with {n{n — l)/2)^ components. This 
representation is reducible, and can be decomposed as 
[2] (8) [2]^ = [n] 0 [n - 1,1] 0 [n - 2, 2]. 

The first term is the invariant = ^^i^Taabb where 
all the indices are contracted. Defining TrQ = Qaa “ 
not to be confused with the symbol tr which corresponds 
to the trace in the fundamental representation - this op¬ 
erator can be expressed as T^(r^) = — 2Tr(5(i)Tr(5(i') 0 
2Tr(3(i)(5(i') = 2{TiQ{i)Q{i') — 1), which is clearly in¬ 
variant. This operator corresponds to a lattice version of 
the energy operator. More precisely, we define the energy 
operator e{ri) = T^{ri) — (T^), where we have subtracted 
the constant term (T^) = 2{n — l)P^^^(i, i') 0 2(l/n — 
l)PQ^^(i, i'), where PQ^^(i,'i') is the probability that i and 
i' belong to different loops, i.e. Po^^(i,'i') = 1 — (i, i'). 

Note that since the sites i and i' lie within a small neigh¬ 
borhood around the position r^, (T^) is indeed a constant 
in a coarse grained picture, z.e., it does not depend on 
Vi. The two point function {e{ri)e{rj)) can be expressed 
exactly in terms of lattice probabilities conditioning the 
points i, i', j and j' to belong to various loops. The ex¬ 
plicit result is not very illuminating, but it shows that 
{e{ri)e{rj)) vanishes as n ^ 1. Denoting by As(n) the 
scaling dimension of £(ri), we thus expect, in the contin¬ 
uum limit 

(£(ri)£(rj)} = Ao(n)(n - ( 7 ) 

where Ao(n) is a regular function of n with Ao(l) ^ 0 
finite, and r = \ri — rj\. In two dimensions, we hav^^ 
A^{n) = I — 1, where g was defined above. 

The second term [n — 1,1] in the decomposition of 
[2] 0 [2]* is the adjoint representation, corresponding 
to the n2 - 1 fields T^.in) = Ec^abcc - 0^° = 
-Q{i)TrQ{i') + Q{i)Q{i') A i ^ i' - This op- 

erator has the same symmetry as (pabi'^i)^ we thus expect 
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its two-point function to be dominated by a power-law 
behavior with the same scaling dimension A 0 (n), with 
some eventual subleading contributions. This correlator 
can also be expressed on the lattice as 


(ri(rOr,yr,)) = -{6ai,-- 


1 \ (n - 2 ) 




+(n - 2){n¥^j\ri,rj) - + ^¥\1!{ri,rj)) 




( 8 ) 


In this expression, ¥2 {^i^rj) is the probability that i 
and j belong to the same loop while i' and j' belong to 
another loop, or that i and j' belong to the same loop 
while i' and j belong to another loop. ¥\ is the 

probability that the four points i, i', j and j' belong to 
the same loop, whereas ¥\/{ri,rj) counts configurations 
in which three out of the four points belong to the same 
loop, while the last remaining point in its own loop (see 
Fig.^. Finally, corresponds to configurations 

where either i or i' is in the same loop than either j or j', 
with the remaining two points being alone in their own 
loop. All these probabilities vanish as r = —rj \ 00 . 

Finally, the last term [n — 2 , 2 ] in [ 2 ] (g) [ 2 ]* has 
dimension n^{n + l)(n — 3)/4, which becomes for¬ 
mally — 1 when n = 1 . As we will describe in 
details in the following, this is the representation 
that will be mixed with the energy operator when 
n = 1. These operators are given explicitly by 

"^abcd ~ "^abcd — {^abT^d + ^cdT^b ~ ~ ^adT^b) ~ 

n(n-i) i^ab^cd “ ^ad^cb) , and they satisfy the con¬ 
straints T^acd ~ point function of then 

reads 


• {•) 


|0 m) 




Dp) 




s(2) 


r 




c5f) 

c;) 



pi,i 

pl,0 

p0,0 


FIG. 3: Graphical representations of loop configurations con¬ 
tributing to the various probabilities introduced in the text. 
Note that some of the probabilities contain more terms than 
the ones drawn: for instance, counts configurations in 
which any of the four points is isolated in a loop while the 
other three points belong to another loop. See main text for 
a precise definition of all these probabilities. 


the scaling limit, we expect 

{TlMTi,,ird) = + 5adShc - ^ x 

i^ac + ^bd + ^ad + ^bc) + 7 - 

(n-l)(n- 2 )y 

( 10 ) 


{'^aabbi'''i)'^ccddi'^j)) ~ ^2 (^ac^bd + ^ad^bc ^ _ 2 ^ 
{^ac ¥ Sbd + Sad ¥ She) + 

(9) 

with a 7 ^ 6 and c ^ d - otherwise this correlator is sim¬ 
ply zero. Despite this rather complicated expression, 
this correlator is remarkably simple. The combination 
of Kronecker deltas is entirely fixed by the symmetry 
constraints ^a=c d'^aacd ~ dependence on 

^2 provides a natural physical interpretation of 

as a “4-leg watermelon operator”, with two loops 
propagating from to Vj. One can also check that all 
the crossed two-point correlation functions between T^, 
and 5 vanish, as required by symmetry. The scaling 
dimension of readJ^ (n) = ^ and in 


3. Other observables 

It is obviously possible to generalize the above to 
observables acting on A > 2 sites, but this will not 
be necessary for our purposes. Moreover, it would 
be natural to ask what happens to the symmetric ob¬ 
servables To^jjQdij^i) Qabid^Qcdid ) Qadi'^^Qebi'^ ) T 
Qab{i')Qcd{i) + Qad{'i')Qcb{i), which would also trans¬ 
form reducibly under SU{n). We chose to focus on 
the antisymmetric tensors Tabcdi'^i) = —Qab{i)Qcd{^') + 
Qadid^Qcbid ) Qabid edid^adid ^Qcbid) to forbid the 
cases a = c OT b = d, which would effectively reduce to 
operators acting on A = 1 site in a coarse grained pic¬ 
ture. It is of course possible to carry on the same analysis 
for Tabcd instead of Tabcdi and one ends up with three ir¬ 
reducible representations with dimensions 1 , — 1 and 

n^(n + 3)(n — l)/4 - note that the dimension of this last 
representation becomes formally 0 at n = 1 (instead of 
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— 1 in the antisymmetric case above), indicating the ab¬ 
sence of mixing. One finds the same pole structure as 
for A/" = 1, and in particular, the limit n ^ 1 of all the 
irreducible operators constructed out of Tabcd is well de¬ 
fined. In the end, we find that the correlation functions 
of these other operators in the limit n ^ 1 do not contain 
any new information so that it is possible to focus only 
on the observables considered in the previous paragraph. 


C. Logarithmic correlations in the n ^ 1 limit 


The n ^ 1 of eq. (10) is clearly ill-defined. This indi¬ 
cates a “mixing” between the scaling operators and 5 
at n = 1, and is further confirmed by the fact that the 
scaling dimensions = At 2 = | coincide for that value 
of n. In order to obtain finite correlation functions for 
n ^ 1, we follovl^ (see alscP^ and introduce a new field 
= (1 - ^ab){T^abb + n(n-i) ^)- Finitcness of correla¬ 
tion functions at n = 1 then requires Ao(l) = ^ 2 ( 1 ), and 
one ends up with the logarithmic correlation function at 
n = 1 


{'ipab{ri)tpcd{rj)) = 2Ao{l)r {SacSbd + Sadhc 

+^ac + ^bd + ^ad + ^ log : (H) 

aJ 

where a is the lattice spacing (UV cutoff), and n is uni¬ 
versal and given by 


/^ = 4 lim 

n^l 


At-2 — Ag 

n — 1 


8^3 

TT 


( 12 ) 


This proves the emergence of a field with logarith¬ 
mic correlations at the critical point. Moreover, 
one can readily check that 'ipab is mixed with the 
energy operator under a scale transformation r 1 -^ 
Ar, a feature characteristic of logarithmic opera- 
torJISl '0^^(Ar) = {i^abir) + f log A 6(r)), whereas 

e(Ar) = A“^/^e(r). We also note that the Kronecker 
delta in 0 are purely formal since n = 1, but it will 
turn out to be convenient to keep track of this structure 
to interpret this correlation function geometrically. 

In order to elucidate the geometrical meaning of this 
logarithmic two-point function, it is fruitful to go back 
to lattice correlation functions. Using eq. (§ and the 
exact expression of the two-point function of 5 in terms 
of probabilities, one can readily express the two-point 
function of ^jjab in terms of various probabilities as well. 
The now well-defined n ^ 1 limit then yields 


{'>Pab{n)'lpcd{rj)) = {SacSbd + ^ad^bc 

-^^ac + ^bd + ^ad + ^bc) + 4 {Fd{r ) — + ^c(^)) • 

(13) 

In this expression, F(i{ri,rj) = — P^’^ is a 

linear combination of disconnected probabilities (Fig. [^: 
P^’^ is the probability that the four points i, i', j and 


j' belong to four different loops, P^4 ig the probability 
that i and %' belong to the same loop while j and / 
belong to another loop, P^’^ is the probability that j 
and 2 ' belong to the same loop while % and %' belong to 
two other loops, or the same with the role of and 
reversed. F(i{ri^rj) has a non-vanishing limit when r = 

|u ^00 given by F^ = (P^^^(i,'i') — P^^^(i, i'))^ = 
(1 — 2P^^^(i, i')Y^ which is a constant independent of the 
position Ti - it is however non-universal and depends on 
how the coarse graining procedure is defined. Finally, 
Fciri^Vj) corresponds to the connected part and is given 

by Fc{ri,rj) = - ¥f,\ri,rj) + ¥^y}{ri,rj) - 

Comparing eq 0 and eq 0, we immediately find 
2¥‘i\r) ~ Ao(l)r-5/2 and 2{Fd{r) - F^ + Fc(r)) ~ 
Ao(l)r“®/^Klog In particular, the ratio 


p(2)(^) a 


(14) 


should be purely logarithmic for a dense loop model with 
fugacity n = 1, with the universal amplitude k given 
by (12). We expect this expression to be particularly 


suited for numerical checks, since it isolates the logarith¬ 
mic dependence. It is worth pointing out that the log¬ 
arithmic dependence stems only from the disconnected 
part Fd{r): using the n ^ 1 limit of (|^, one can easily 
show that the connected piece Fc{r) oc (T^^(r)T^\(0)) is 
purely algebraic for n ^ 1 , with a leading contribution 
given by the 2 -leg exponent A 0 (n = 1 ) = |. 


D. Relation to the Potts model, percolation, and 
supersymmetry 


1. Percolation and Potts model 

We have thus shown that dense loop models with fu¬ 
gacity n = 1 involve a logarithmic operator resulting from 
the mixing of the energy and the 4-leg watermelon fields. 
The two point function of this logarithmic operator can 
be expressed in terms of probabilities that could readily 
be measured using Monte Carlo simulations for example. 
When n = 1, the two dimensional dense loop model in¬ 
troduced above can be mapped onto bond percolation, 
where the loops are the hulls of percolation clusters. Our 
results thus also apply to the two-dimensional percola¬ 
tion problem, where the observables that we described 
act on the hulls of percolation clusters. The existence of 
a logarithmic operator mixing energy and 4-leg fields in 
percolation was pointed out in Refs. I28|5Q| and a similar 
concre te log arithmic observable in percolation was un- 
coverecP^®! using the Q ^ 1 limit of the Potts model. 
Despite the similarities between the above calculation 
and the result of Ref. [28] - that can be traced back to 
the relation between the representation theory of SU (n) 












and the symmetric group - note that the resulting loga¬ 
rithmic observables are different, since the ones obtained 
from the Potts model act on percolation clusters, instead 
of acting on hulls. Therefore, even if the Potts model is 
known to be related to loop models with SU(n) symme- 
tr}!^, the observables in both models are different. For 
instance, the probability that two points separated by 
a distance r belong to the same percolation cluster is 
definitely different from the probability that two points 
belong to the same loop surrounding a cluster, and they 
actually sc ale wi th different exponents, and 

respectivel}@^^. 


Supersymmetry vs Replieas 


Although the loop models introduced above do not in¬ 
volve quenched disorder, it is helpful to think of our ap¬ 
proach based on the formal limit n ^ 1 of the SU{n) 
symmetry as a sort of replica trick. As in disordered 
systems, it is possible to use an alternative, more rigor¬ 
ous ap proach based on supersymmetry to construct loop 
modelP^®^. Intuitively, one introduces both fermonic 
and bosonic operators to ensure that closed loops get a 
weight n = 1, the fermonic terms giving negative contri¬ 
butions when computing the weight. Loop models with 
fugacity n are then described by a SU{n + m|m) global 
(super)symmetry (SUSY), where SU{n-\-m\m) is the su¬ 
pergroup analog to SU{n)^ defined by the transforma¬ 
tions preserving a form with n^m bosonic variables and 
m fermonic ones (see e.g. Ref. [53] for a review of Lie su¬ 
peralgebras) . There is then no iss ue wi th the value n = 1 
and percolation can be describecP^^^ by, e.g. a theory 
with SU{2\1) SUSY (or SU{m^l\m) SUSY in general). 
The emergence of logarithmic correlations can also be 
understood in this supersymmetric language !^ us¬ 
ing more involved representation theory tools, such as in¬ 
decomposable representations - representations that are 
reducible but not fully reducible. We will not describe 
the SUSY approach in more detail here (see EQl), but we 
simply note that although SUSY is very powerful to pre¬ 
dict formally the emergence of logarithmic operators, it 
does not yield simple geometrical formulas like (14). 


3. Extended symmetry 


As mentioned above, the actual symmetry at the 
strong coupling fixed point of the two-dimensional 
sigma model is much larger than SU{n^^. As a result, 
operators that transform under different irreducible rep¬ 
resentations of SU (n) could actually be part of the same 
multiplet of fields for the actual, larger symmetry. This 
means that the prefactors of some of the logarithms pre¬ 
dicted using the SU(n) symmetry could in fact be zero. 
In the example 0 described above, the prefactor 
is non-zero and can be computed exactly. Moreover, the 
true symmetry of the critical point is also known exactly 


for those dense loop modelP^. For example, the repre¬ 
sentation [n — 2,2] with dimension rU{n + l)(n — 3)/4 
of SU{n) is part of a much larger representation of di¬ 
mension rU{nS — 3)/2 for the actual symmetr^ESl (or 
rU — 3rU + 1 for fields living at the boundary). How¬ 
ever, in the limit n ^ 1, the dimension of both repre¬ 
sentations goes (formally) to —1, indicating a “mixing” 
(in the sense described above) with the trivial represen¬ 
tation of dimension 1. Therefore, in the limit n ^ 1, the 
SU(n) symmetry is enough to understand the mixing of 
the 4-leg operator (corresponding to the representation 
[n — 2, 2]) and the identity operator. 

Even if in this case, considering a smaller symmetry 
does not affect the prediction of logarithmic correlations 
for the 4-leg operator, this example nicely illustrates that 
fields having different symmetry properties under SU(n) 
can actually have the same scaling dimension. More pre¬ 
cisely, for bulk fields we find that the large represen¬ 
tation n^(n^ — 3)/2 for the true symmetry decomposes 
as (n^(n + 3)(n — l)/4) 0 (n^(n — 3)(n 0 l)/4) under 
SU{n)^ where we denoted representations by their di¬ 
mension for simplicity. For boundary fields, the repre¬ 
sentation rU — 3n^ 01 for the true symmetry decomposes 
as (n^(n03)(n-l)/4)0(n^(n-3)(n0l)/4)0(n^-l)02x 
((n^ — 4)(n^ — l)/4). These irreducible SU{n) represen¬ 
tations - all contained in the tensor product ([1] 0 [1]*)*^^ 
- therefore correspond to the same primary field (the 4- 
leg operator) for the extended symmetry. Interestingly, 
some aspects of that analysis can be recovered directly 
at the level of the sigma modeP^. 

Note also that considering a smaller symmetry group 
at the critical point may also lead us to miss some log¬ 
arithms: for example, considering the full symmetry of 
dense loop models predicts a mixing between the repre¬ 
sentations of dimension n^(n^ ~ 3)/2 (or — 3n^ 0 1 at 
the boundary) and rU — 1 (adjoint) for n = \/2, consistent 
with other approache^^. However, we would completely 
miss this mixing by considering only the smaller SU{n) 
symmetry - there is no pole at n = y/2 in eq. ( p^ . Hav¬ 
ing these drawbacks in mind, we are now ready to go 
back to Anderson transitions. 


III. LOGARITHMIC OBSERVABLES AT THE 
SPIN QUANTUM HALL TRANSITION 

We now turn to the study of logarithmic correlations 
at the Spin Quantum Hall transition, a close cousin of 
the Integer Quantum Hall transition where SU (2) spin is 
transported instead of U{1) charge. We first review the 
basic physics of the Spin Quantum Hall Effect and in¬ 
troduce the corresponding network model that describes 
the transition. Using the replica trick and the global 
symmetries at the critical point, we then exhibit linear 
combinations of Green’s functions whose disorder average 
should be logarithmic. Einally, using supersymmetry and 
an exact mapping onto classical percolatiorP^, we argue 
that this logarithmic observable found using the replica 
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trick coincides with the logarithmic two-point correlator 
in the SU{n = 1) loop model discussed in the previous 
section. 


A. Spin Quantum Hall Transition 

The Spin Quantum Hall Effect (SQHE) is the quanti¬ 
zation of the spin Hall conductance, analogously to the 
quantization of the Hall charge conductance in the Inte¬ 
ger Quantum Hall Effect (IQHE). Spin transport is de¬ 
fined as the response to the spatial variation of an exter¬ 
nal Zeeman magnetic field, which plays the role of the 
electric field in the IQHE. If the external field is along 
the 2 : direction, the spin Hall conductance a^y is defined 
by the spin current along the x direction 

Jx = (15) 

As in the IQHE, the spin Hall conductance a^y can take 
non-zero values only in systems with broken time-reversal 

invariance. _ 

The SQHE may occur is singlet superconductor^^^^ 
with mean-field Hamiltonian 

ij a ij 


There is no notion of charge transport since charge is 
obviously not conserved by the superconductor, but the 
Hamiltonian has a global SU (2) spin rotation symmetry 
if Aij = Aji. If the gap function A is complex, time 
reversal invariance is broken and a^y may take non-zero 
values. Eor a clean dx 2 -y 2 ^id^y superconductor, one can 
show that (T^y is quantized and that transport is carried 
by chiral edge modes at the edge of a system with a 
boundary. These edge states are robust to weak disorder, 
whereas strong disorder eventually drives the system to 
a topologically distinct localized phase characterized by 
(j^y = 0 (spin insulator). In the following, we will study 
this Anderson transition from the SQHE phase to the 
spin insulating phase. 

Spin rotation symmetry and broken time reversal in¬ 
variance defines the symmetry class Th e corr e- 

sponding sigma model describing the transit iorl^^^^ is 
defined on the coset Sp{2n)/U{n) with n ^ 0 using a 
replica approach, or OSp{2\2)/U{l\l) using supersym¬ 
metry (more general cosets OSp{2m\2m)/U{m\m) can 
also be used to describe fluctuations of observables and 
compute multi-fractal exponents). 


B. Network model and Green’s functions 

The lattice version of this non-interacting di sorder ed 
fermions problem is introduced using a varianfP^^ of 
the Chalker-Coddington network model for the Integer 
Quantum Hall transitiorP^. 



FIG. 4: Network model with directed links and two types of 
nodes A and B on the square lattice. Particles can propagate 
along the directed links, with random SU{2) (resp. U{1)) 
matrices associated with the propagation on each link in the 
SQHE (resp. IQHE) case. The blue and red circles represent 
scattering matrices at each node. 


As in the Chalker-Coddington model, we introduce a 
network with directed links and two types of nodes A and 
B on the square lattice (Eig. ^ where particles can prop¬ 
agate. The difference with the Chalker-Coddington case 
comes from the spin of these particles, as disorder is in¬ 
troduced by random SU (2) matrices (instead of random 
1/(1) phases) associated with the propagation along each 
link. The nodes of the network correspond to scattering 
events and we take the scattering matrices at each node 
to be diagonal in the spin subspace and 


SK,a 


y 1 - t'jia tKa \ 
~tKa yi ~ t'^Ka ) 


(17) 


where cr labels the spin of the particles and K = 

A^B denotes the two different types of nodes. Spin 
rotation symmetry implies that Ik = = ^^ 4,7 

and isotropy of the network enforces 't\ ^ = 1 , so 

that we parametrize Ia = cos 6 > and Ib = sinO. The 
isotropic critical point of the network model then occurs 
for = ts = and varying 0 away from 7 r /4 drives 
the system into either a spin insulator or a quantum spin 
Hall state, with a jump in the spin quantum Hall con¬ 
ductance at the transition. 

The simplest way to define Green’s functions is to use 
a first quantized formalism where the single particle wave 
function ' 0 (r, t) follows a discrete time evolution ' 0 (r, t + 
1 ) = U'ilj{rU) where the evolution operator // is a sparse 
2N X 2N symplectic matrix for a network with N links. 
U simply describes the discrete evolution of the wave 
function with random SU{2) matrices on each link and 
the correct scattering factor at each node. The retarded 
Green’s function between two links i and j is then defined 
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as 


zU) (18) 


where 2 ; = with the energy e = 0 at the transi¬ 

tion and 77 = 0+ an infinitesimal level broadening. The 
advanced Green’s function is defined by replacing ^ by 
but because the evolution operator is a symplec- 
tic matrix, it is possible to relate the advanced Green’s 
functions in terms of the retarded one, so that all phys¬ 
ical observables can be expressed in terms of the re¬ 
tarded Green’s function (18) only. At the transitiorP^, 


the Gree n’s functions sh ow a critical behavior, so that for 
example trG(i, j)G(j, i) ^ Note that the averaged 

local density of states is critical and scales as p{e) ^ 
at the transitioiP^, in sharp contrast with the IQHE case. 


C. Second quantization and supersymmetry 

The SQHE network model can also be conveniently de¬ 
scribed within a second-quantized framework. One then 
thinks of the evolution operator as a product of trans¬ 
fer matrices, and introduces fermions to represent the 
retarded Green’s function. Usually, one needs both re¬ 
tarded and advanced operators to describe the product 
of retarded and advanced observables, but as already de¬ 
scribed above, we will only need the retarded Green’s 
function in the case of the SQHE effect. In order to 
average over disorder, we use either the replica or the su¬ 
persymmetry trick. Within the replica trick, the Green’s 
function can be represented as 

= nm(/“(z)/“t(j)), (19) 

with, in the path integral formalism, 

(/“(*)r|(i))= / 

( 20 ) 

where = /a (^) ^ vector of size 2N, 

and the path integral measure reads P[/, /] = 

na=i rii cTi Because of the 

replica limit n ^ 0 , the average over disorder can be 
realized explicitly and it is useful to think of the critical 
point of the disordered system (the SQHE plateau transi¬ 
tion) as obtained from a limit n ^ 0 of a family of GETs 
with central charge c{n) ^ 0 as n ^ 0 . 

A more productive approach in the context of the 
SQHE is to use the so-called supersymmetry trick in or¬ 
der to average over disordeil^ (see also Refs. l6QI6ip . In¬ 
stead of a family of fermions (or bosons) parametrized by 
n ^ 0 as in the replica tick, one introduces two bosons ha- 
and two fermions fa per link (one per spin), with canon¬ 
ical commutation relations on “up links” of the network, 
and modified anti-commutation relations for the fermions 
on down links (see Ref. [36] for details). The equal num¬ 
ber of fermions and bosons ensures that closed loops in 


the graphical expansion of the partition function van- 
ish, so that ^susY = 1- In the path integral formal- 
isnpil, this can be understood very easily as ^susy = 
/T>[/,/]P[ 6 ,^ thanks to the Gaus¬ 


sian integral identities f P[/,= Det(l — zU) = 
(f Vlb, In the second quantized framework, 

one ends up with a transfer matrix that commutes with 
the superalgebra osp(2|2) (see e.g. Ref.|53jfor a review of 
Lie superalgebras) for any disorder realization, and the 
average over disorder then projects the Eock spaces on 
each up or down links onto SU(2) singlets. The resulting 
Eock spaces can be shown to correspond to the fundamen¬ 
tal representation of sl( 2 |l) ^ osp( 2 | 2 ) (for up links) and 
its dual (for down links), both with finite dimension 3 (su¬ 
per dimension 2 | 1 ) with two bosonic states and one fer- 
monic one. This drastic simplification - the Eock spaces 
become finite-dimensional after averaging over disorder 
- allows one to solve this supersymmetric model exactly 
by mapping it onto the classical percolation problem. We 
will not need the details of this mapping here and instead 
refer the interested reader to Ref. EH but we simply re¬ 
mark that using the isomorphism osp( 2 | 2 ) sl( 2 |l), the 

resulting loop models obtained by expanding graphically 
the partition function corresponds precisely to the super- 
symmetric “5'U(2|1) ^ SU{ n ^ 1 ) ^ percolation” model 
studied in Sec. [H| (see [IIP 2 in particular). This mapping 
is summarized in Eig. Note that the only SQHE ob¬ 
servables that can be mapped onto classical percolation 
in this way can be described using two bosons ba and 
two fermions fa per link only. Eor instance, multi-fractal 
properties usually require more species of fermions and 
bosons, with supersymmetry osp( 2 m| 2 m) and m > 1 , 
and this mapping onto percolation breaks dowrl^^^. 


D. Logarithmic correlations from the replica trick 


1. Observables aeting on N = 1 link 


We now focus on the replica approach and start from 
fermions fa = fa^ with a = (a, cr) and a = l,...,n, 
with n the number of replicas. We work with fermionic 
operators because bosonic variables in Anderson tran¬ 
sitions typically involve non-compact symmetries, lead¬ 
ing to a much more complicated representation theory 
analysis. Eor symmetry class G, we expect the global 
symmetry of the network model to be 5'p(2n), the group 
of unitary matrices U satisfying = cr^. We start 

with observables of the form /a(i)/^(i) where i labels the 
links of the networks. The conjugate of the fundamen¬ 
tal representation [ 1 ]* of Sp{2n) is pseudo-real [I]"*" [ 1 ] 

since U* = a^Ua^. It is therefore convenient to intro¬ 
duce wMch transforms as c 1 -^ Uc 

if f Uf with U G Sp{2n). The combination f'^cr'^c is 
then clearly invariant, and the n(2n — 1 ) — 1 antisymmet- 
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FIG. 5: Percolation description of the SQHE transitiorP^, with periodic (resp. open) boundary conditions in the vertical 
(reps, horizontal) direction. The blue loops can be thought of as (semiclassical) paths contributing to spin transport, and 
their orientations correspond to alternating fundamental and dual representations of the superalgebra sl(2|l). Average over 
disorder configurations enforces that these loops be non-intersecting, so that they can be mapped onto hulls of classical 
percolation clusters (shown in red in the middle picture only). The SQHE transition then corresponds to the critical point of 
the percolation problem, separating a trivial spin insulator phase from a SQH insulator phase with a jump is the spin quantum 
Hall conductance. 


ric operators 

y 

Wa/3 = faC/3 - f^Ca + ( 21 ) 


which is ill-defined in eq. ( [2^ . The only way to fix this 
divergence is to require A( 0 ) = 5(0) and A 0 (O) = A<^( 0 ), 
and introduce = 0 aa + The n ^ 0 limit of is 
then well-defined 


form an irreducible representation of Sp{2n). The n( 2 n+ 
1 ) symmetric counterparts V^a /3 = /aC /3 + f^Ca then cor¬ 
respond to the adjoint representation of Sp{2n). The 
average over disorder in the network model enforces the 
physical operators to be SU ( 2 ) singlets. Among all these 
operators, we therefore restrict the physical observables 
to 


{'ipa{ri)'ipbirj)) = A(0)r 2 A<^(o) + ^log ^) , (24) 

with K = 21 im^^o(A 0 — A$)/n. The same correla¬ 
tion function can also be expressed in terms of disorder- 
averaged Green’s functions in the SQHE network prob¬ 
lem using Wick’s theorem 


Mi) = = E /-“/r - V E ’ (22) 


C=1 


which transform under the adjoint r epresen tation of the 
subgroup U{n) C Sp{2n) (see Sec. HID 2 | below). We 
thus expect the two-point function of (pab for generic n 
to take the form 

{Mn)Mrj)) = Ain) (^ 5 ab - ( 23 ) 


where the scaling dimension A^^ is obviously differ¬ 
ent from the one used in Sec. |l^ for SU{n) models, 
in particular, its general expression as a function of n 
is not known. We also define the invariant operator 
^ = Ea Ec=ialong with l> = $ -($). The 
two-point function of is then given by {^{ri)^{rj)) = 
with 5(0) ^ 0 and A<^(n) is the scaling 

dimension of 

Disorder-averaged correlation functions in the network 
model are then obtained by considering the limit n ^ 0 , 


{'Ipa{n)'ipbii'j)) = tiGii,i)tTG{j,j)-trGii,i)xtTG{j,j) 


- SabtTGii,j)Gij,i), (25) 

where tr represents the trace over the spin index - recall 
that the retarded Green’s function G is actually a 2 x 2 
matrix. Comparing equations (24) and (25), we can then 
infer the scaling of various Green’s functions as a function 
of the distance r. 

As discussed above, using the supersymmetry map- 


ping^, one can then relate these Green’s functions to per- 

. For instance. 


eolation observables (see also Refs, 
it is straightforward to show that trG(i, i)trG(j, j) = 
E.„.,(^a.(i) 6 t,(i)/.,(j)/J,(j))susY = ((1 + 2 B(i))(l + 
SVSY = 1, where 5 and Qz are two of the eight 
SUSY generators of sl( 2 |l) in the fundamental represen¬ 
tation (see e.g. Ref. [53] or the appendix of Ref. [63|) and 
we have used the fact that sirBQz = str5 = strQ;^ = 0, 
with str the su pertrace in the fu n damenta l represent a- 
tion. Therefore, trG(i,'i)trG(j, j) — trG(i, i) xtrG(j, j) = 
0 so the correlation function that could show a loga¬ 
rithm at this order is trivially zero. Similarly, one finds 
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that the non-logarithmic part trG(i, j)G(j, i) is given by 
—where is the probability - intro¬ 

duced above in Sec. [Ill “ t hat i and j be long to the 
same loop. This shows^that trG(i, j)G(j, i) decays with 
the 2 -leg ( 1 -hull) exponent which allows us to identify 
A 0 (n = 0) = J. These results imply that the param¬ 
eter K in eq. ( [^ is actually zero. Note that nothing 
in our symmetry analysis prevents this from happening. 
This could be just a coincidence, i.e. the derivative of 
the exponents A 0 (n) and A<|,(n) with respect to n could 
coincide at n = 0 , but this could also indicate the pres¬ 
ence of a larger symmetry that enforces A0(n) = A^{n). 
We emphasize that the actual symmetry of the critical 
point of the SQHE is most likely must larger that U(n) 
or Sp{2n), so we always run into the risk of predicting 
logarithms that appear with amplitudes /i: = 0. This was 
also true for the loop models of Sec. [TT|- for which it is 
known that the actual synimetry of the critical point is 
much larger than SU{n]^^ but there we had a way to 
compute K exactly, see eq. ( 12 ). Actually, had we consid¬ 
ered the above SU(n) model at a mean-field level (above 
the critical dimension), we would have also found /^ = 0 
instead of ( [T^ . Our approach is nevertheless fully con- 
sistent - the s ymmetry analysis correctly predicts that 
trG(i, j)G(j, i) should be a scaling operator - and we de¬ 
duce that observables acting on a single link do not have 
logarithmic correlations. 


Sp{2n) vs U{n) symmetry 

In order to find logarithmic observables, we thus turn 
to operators acting on two nearest neighbor links i and i' 
as in Fig. Before doing so, we first remark that equa¬ 
tion ( [^ suggests restricting to S'!/( 2 )-invariants oper¬ 
ators from the very beginning: analyz¬ 

ing the U{n) symmetry of the remaining indices. Inci¬ 
dentally, the remaining symmetry U{n) coincides with 
the denominator in the coset Sp{2n)/U{n) of the sigma 
model describing the transition. This simplification is 
motivated by the fact that the genuine symmetry of the 
critical point is unknown anyway, so that we choose to 
work with the subgroup U(n) of Sp(2n) in order to sim¬ 
plify dramatically the calculations. As we will see in 
the following, this will turn out to be enough to iden¬ 
tify observables that should have logarithmic correlation 
functions at the critical point. 

We remark that this is related to our choice of con¬ 
sidering quadratic observables /a(^)/^(^) from the begin¬ 
ning. It is of course quite natural starting from operators 
like /a(^)/^(^) that the SU{2) invariance has to be en¬ 
forced within the replica approach before taking the limit 
n ^ 0 since the spin structure is completely lost in the 
replica limit Sp{2n ^ 0 ). A more correct approach would 
be to consider more complicated representations V and 
of Sp{2n) acting on a single edge, that are compati¬ 
ble with the SU ( 2 ) symmetry and that involve products 
of fermionic operators with different replica indices (see 


Ref. [38] in the context of the IQHE). Other observables 
would then be constructed by taking tensor products of 
these representations. However, we will see in the follow¬ 
ing that some of the logarithmic structure of the SQHE 
can be understood by starting with the much simpler 
quadratic operators fafV satisfying the SU{2) local 
(“gauge”) invariance, classifying them using the remain¬ 
ing unitary symmetry of the replica indices and analyz¬ 
ing the limit n ^ 0. We also note that in the SUSY 
language, this amounts to considering the much simpler 
subgroup U(I|I) of OSp{2\2) ^ U( 2 |I) whose indecom¬ 
posable representations are suspected to play a crucial 
role in disordered systemJI^. 


3. Observables aeting on N — 2 links 


We now turn to observables acting on two nearest 
neighbor sites i and i'. Let us consider operators of the 
type which are mani- 

festly SU{2) invariant. As in Section [Ilj we expect ob¬ 
servables symmetric under the exchange a ^ c or b ^ d 
to reproduce the physics of the operators acting on A = I 
link, with no logarithm for n = 0. We thus focus on 
the antisymmetric tensor Tahcd{ri) = -Qah{i)Qcd{i') + 
Qad{d)Qcb{d ) Qabid ^Qcdid) Qadid ^Qcb{d)'i with 
Qabii) = Wc then cousidcr = 

^ab'^aabb which is invariant under I/(n), and introduce 

$( 2 ) ^ ^( 2 ) _ ^$( 2 )^_ rpjjg £gljjg ^^ 2 ) ^ Y.c'^abcc then 

transforms in the adjoint representation of U{n). These 
two operators and have the same symmetry as $ 
and (j)ab introduced in sect ion [TlIP I[ and the same mech¬ 
anism leading to logarithms in the limit n ^ 0 applies 
here as well. Because they have the same symmetry, the 
ill-defined limit n ^ 0 indicates logarithms in subleading 
contribution^^, so that we find, in the limit n ^ 0 


= A8abr 

+ (5^6 + C log ^) + ..., (26) 


where ip^'^ = piaa + we have seen above, 

^ 0 ( 0 ) = \ corresponds to the one-hull percolation expo- 

nent, while A^ ^(0) is a priori unknown from this replica 
analysis, but we will shortly see that it is given by the 
two-hull percolation exponent, A^^^(O) = |. 

To translate this into a concrete prediction for the net- 

(2) (2) 

work model, we evaluate the correlator {^pa vu )'05 A^i)) 
in terms of disorder-averaged Green’s functions. We find 
16 different contributions which we compute using Wick’s 
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theorem, including for example 


n 

]im 5I(/a,<Ti(Pl)/cVi(^^l)/c,<T2(P2)/i,<,2(^^2) 

c,d=l {cr} 

X A,<73 (P3)/],<,3 (P3)/d,<74 {Pi)flai (P4)) 

n _ 

= i™0 51 51 i 4 (P* ’ Pi ) > (27) 

c,ci=l{a} ’■■■’ 


with (a G {a,c, 6, d}, (3 G {c, a,(i, 6} and p G j, J^}- 
Computing these 16 determinants explicitly and taking 
the limit n ^ 0, we find an expression that can be recast 
as 


{i’a{ri)tpb{rj)) = 4 (n(r) - Ooo + 5(r) + ©(r)) 

- 46ab (2 X S(r) + 0(r)). (28) 


The most important piece of this correlator is the discon¬ 
nected part 

n{r) = tTG{i,i)tvG{i' ,i')tTG{j, j)tvG{j', j') 

+ tr [G{i, i')G{i', i)] trG(j, j)trG{j',j') 

+ trG(i,i)trG(i',P)tr [G{j,j')G{j',j)] 

+ tr \G(i,i')G{i', z)] tr \G{j,j')G{j',j)], (29) 

which has a finite limit as r ^ oo: 

^oo = ^trG(i,'i)trG(i', i') + trG(i, . (30) 


The functions 0(r) and S(r) can also be readily com¬ 
puted and contain 16 and 4 terms respectively, but their 
explicit expression will not be relevant to our purposes. 
Indeed, in analogy with the discussion in Sec. [ll| we ex¬ 
pect the logarithmic correlations to arise from discon¬ 
nected terms such as 11 (r), and not from connected ob¬ 
servables like 0(r) and S(r) - at this point, this state¬ 
ment is a conjecture based on analogy but it will be ver¬ 
ified below when we compare our results to the super- 
symmetric approach. Based on this and using eqs. (26) 
and ([ 2 ^, we thus obtain 


H(r) — Hoo ^ + /3r log - + 


(31) 


with a the lattice spacing (UV cutoff), where we have 
anticipated the fact that (0) = | (see supersymmet¬ 
ric approach below). Even though is enough to show 
the existence of logarithmic correlations in the disordered 
SQHE network model (which would naively contradict 
scale invariance and therefore immediately implies the 
existence of logarithmic operators!), it may seem like a 
rather weak prediction since the logarithm appears only 
as a correction to a power-law function, which itself is 
a subdominant contribution. We will see in the follow¬ 
ing that the IQHE transition leads to a prediction that 
should be (much) more appropriate for numerical checks. 


E. Logarithmic correlations from the 
supersymmetry trick 


We now argue that the logarithmic correlation © 
obtained from the replica trick can be recovered using 
an independent method based on the results of Sec. [HI 
and the SUSY description of some SQHE observables in 
terms of percolation probabilities^ (see also Refs. [6Tti6^ . 
Recall that the SUSY description involves two bosons ba- 
and one fermion / on each link, corresponding to the 
fundamental representation of the superalgebra sl(2|l). 
In that language, eq. ( [^ can be written as 

^(^) = Y (^<7i^<7i(*)’''^'<T,,&<7,,(i')^^^<7,(j)’''&<7^<^/(i')''')d, 
fcrcr'} 

(32) 

where r = — r^l, and the subscript d in the correlator 

refers to the disconnected part (with a finite limit r 00 ) 
involving Wick contractions between i and i' or j and j' 
only. This correlation function can easily be expressed in 
terms of the SUSY generator B = + 1)? so 

that it can be expanded onto various classical percolation 
probabilities, with amplitudes given by the supertrace 
of powers of B in the fundamental representation (see 
e.g. Ref. [63] for related calculations). Similarly, we find 

that Uoo = ^trG(i, i)trG(i', i') + trG(i,'i')G(i', i)^ = 

(l+ 4 (B(i)B(z'))susY)" = F^, where = (P^(i,z')- 
P 4 °^(i,i'))^ was introduced in Sec.|l^ Gathering these dif- 
ferent pieces, we find the following expression in terms of 
percolation probabilities 


n{r) - Oo, - F^r) - F^, 


(33) 


where Fd{r) is precisely the linear combination of discon¬ 
nected percolation proba bilitie s introduced after eq. (13). 
As we have argued in Sec. II C[ Fd{r) has a leading power- 
law behavior associated with the 1-hull percolation expo¬ 
nent, with logarithmic corrections to the subleading term 

- associated with the 2-hull or 4-leg percolation exponent 

- that can be understood as the singular limit n ^ 1 of 
SU{n) loop models. Using these results, we thus con¬ 
clude that eq. (31) holds, in agreement with the replica 


analysis. Therefore, the logarithmic correlations in per¬ 
colation described in Sec. m have the same physical 
origin as the logarithm in (31), and they correspond to 


the “mixing” of the energy and 2-hull (4-leg) operators 
in the 2D percolation problem. 


IV. (TENTATIVE) GENERALIZATION TO THE 
IQHE TRANSITION 

We have shown above that both the replica trick and 
the supersymmetry trick - combined with the analysis of 
SU{n 1) loop models - predict the existence of loga¬ 
rithmic correlations at the SQHE plateau transition, as¬ 
sociated with concrete disorder-averaged observables (see 
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eq. ([^). Although the SUSY mapping onto percolation 
is obviously very specific to the SQHE transition, our 
analysis of the n ^ 0 limit of the replica trick can be 
generalized to other symmetry classes. To illustrate this 
point, we conclude this paper by shortly describing the 
case of the Integer Quantum Hall Effect (IQHE), with 
symmetry class A. The network model describing the 
transition was introduced in the seminal work of Chalker 
and CoddingtorP^, the only difference with the SQHE 
case being the random U{1) phases - instead of SU{2) 
matrices - on each link. The associated topological sigma 
modePI with bosonic fields reads 

jC = ^tr(a^(3)^ - ^e^^trQd^Qd^Q, (34) 

where Q lives on the coset U{n^n)/U{n) x U(n) (or 
U(l,l|2)/U(l|l) X U(l|l) in the SUSY formulatioiP^ , 
with n the number of replicas, and axy plays the role of 
a topological angle associated with the non trivial sec¬ 
ond homotopy group 7 r 2 = Z of the target space. At 
the plateau transition, the sigma model flows to strong 
coupling and we expect the quantum critical point to 
have a large symmetry, and we will assume that U(n, n) 
is a subgroup of this symmetry. To avoid having to 
deal with non-compact symmetries, we restrict ourselves 
to the case of fermonic replicas, with symmetry group 
U(2n). We thus write the retarded and advanced Green’s 
functions as G^{i,j) = lini„^o(/^where a = 
1 ,..., n, with the U ( 2 n) symmetry acting on the indices 
OL = (a, zb). 


A. Observables acting on A = 1 link 

Using the symmetry analysis of Sec. |TI| we intro- 

duce the operators (paji = - Y S 7 

which could potentially lead to logarithms in the limit 
n ^ 0 because of the second term. However, the average 
over disorder in the Chalker-Coddington model leads to 
observables with the same number of retarded and ad¬ 
vanced particles on a given link - this can be understood 
in a graphical expansion as retarded paths come with 
a random phase e*^ while retarded paths are weighted 
with integrating over 7 then restricts the number of 
paths of each types to be the same. We therefore consider 
the observables (pabi^) = which transform 

irreducibly under the reduced symmetry U{n) x U{n). 
Note that this is similar to the reasoning that led us to 
consider a reduced U{n) symmetry for the SQHE tran¬ 
sition in order to ensure that all observables were SU ( 2 ) 
singlets, a condition enforced by the average over disor¬ 
der. Once again we emphasize that this is related to our 
choice to start with simple quadratic operators: the ver¬ 
tex model corresponding to the U ( 2 n) /U (n) x U (n) sigma 
model is in fact built out of more complicated represen¬ 
tations V and V'^ of U ( 2 n) (involving products of up to n 
fermion operators) compatible with the local U(l) invari- 
anc^^. As in the SQHE case, we choose to focus instead 


on the much simpler operators 0 a 6 (^) = 
to restrict our analysis to the subgroup U{n) x U{n) of 
U(2n). (Erom the SUSY point of view, this means that 
we are considering only some part of the indecompos- 
ability of the supergroup U(l, 1 | 2 ) corresponding to the 
much simpler subgroup U(l|l)xU(l|l).) We expect this 
analysis to be enough to uncover potentially logarithmic 
observables as n ^ 0 . 

There is no 1/n pole in correlation functions at this 
order (this differs from the SQHE case) and the only 
non-trivial correlator is 

_ (35) 

consistent with the fact that j)G^{j^i) should be 

a scaling operator. 


B. Observables acting on A = 2 links 

As for the SQHE, one needs to consider coarse 
grained observables acting on more than one link 
to generate logarithmic correlations. We thus con¬ 
sider products of operators acting on a single link 
of the form which can be 

thought of as being obtained from ‘fusing’ A = 1 
operators - we recall that i and i' are close neigh¬ 
bors. These operators form a reducible rep¬ 

resentation of U{n) X U{n). There is a unique 
invariant = 

TrQ^ {i)Q ^(^0? where the trace symbol corresponds 
to a sum over the indices and 

There are also two representations corresponding to 
the products trivial x adjoint or adjoint x triv¬ 
ial in retarded/advanced spaces, given by = 

Ec/c'^(*)/a“^(0/r(*')/7(^0 - and = 

- 7"T(0-0). Because of the 

1/n terms in these operators, these operators could po¬ 
tentially have an ill-defined n ^ 0 limit that would even¬ 
tually lead to logarithmic correlations. However, it is 
not hard to see using Wick’s theorem that because of the 
sum in the first terms, the n ^ 0 limit of, say, the cor¬ 
relator j')) is well defined. In other 

words, the correlation function j')) 

will involve a term 5ah — 7 dictated by representation 
theory that, similarly to other examples we have stud¬ 
ied in this paper, could lead to logarithms in the limit 
n ^ 0 ; but it appears with an overall 0{n) pre-factor 
that cancels the 1/n pole. After subtracting the in¬ 
variant and these two representations, we are left with 
— I — 2 X (n^ — I) = (n^ “1)^ operators that trans¬ 
form under the product of adjoint representations under 
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U{n) X U{n). These operators are given by 


Gathering these different pieces, we conclude that 


i’abcd = \i)fc 

^ k 

- (36) 


We thus expect the two-point function 
generic n 7 ^ 0 to scale algebraically as 
an amplitude proportional to {5ad — \){^hc ~ \)' Simi¬ 
larly to the previous examples encountered in this paper, 
we introduce a new operator ^abba = '^abba — (V^) and 
^abba = ft{i)fb^{i)fb^')ft^^') = ^abba + + 

ia order to solve the ill-defined limit 
n ^ 0. Because there are three different operators i^abba^ 
^ T^d Q^] 2 d involved, the resolution of the ‘n ^ 0 
catastrophe’ is more intricate than the other examples 
encountered in this papeJ^. We leave the detailed un¬ 
derstanding of this n ^ 0 limit for future work, but sim¬ 
ply notice that regardless of these details, the n ^ 0 
limit will yield (' 0 a 66 a' 0 cddc) ^ logr if a 7 ^ d and 

b ^ with the possibility of having log^ r terms as well. 
The important point is that in analogy with the other 
examples described in this paper, we expect to have log¬ 
arithmic corrections in the disconnected part of the cor¬ 
relation function ('0a66a'0cddc)’ Obtained by enforcing all 
the replica indices a ^ d and 6 7 ^ c to be different. Using 
Wick’s theorem, we find 

i'lpabba'lplddc) = = h - '^j\) - ^ 00 , (37) 

for a 7 ^ d and 6 7 ^ c with 

n(r) = G+{i,i')G-{i',i)G+{j,j')G-{j',j). (38) 

The function U(r) has a (non-universal) finite limit as 
- 2 

r ^ 00 , given by Uoo = {i 'which is inde¬ 

pendent of Vi - recall that i and i' are chosen to be in the 
infinitesimal neighborhood of r^. The replica approach 
thus predicts 


0(r) — Uo 




-2A^(0) 


log“ r, 


(39) 


with a = 1 or a = 2. Contrary to the SQHE case 
logarithms appear in the leading power-law contribution, 
and there is no need to subtract a very complicated com¬ 
bination of Green’s functions to isolate the logarithmic 
term. The downside is that the limit n ^ 0 is less con¬ 
trolled in the IQHE case as it involves three operators, 
leading to the indetermination of the exponent a = 1 
or (a = 2 in (39). Note also that the term coming with 
amplitude dadhc in the correlation function (' 0 a 66 a' 0 cddc) 
should scale purely algebraically, so that 


r(r) = G+{j,i)G+{i',r)G-{i,3)G-{j',i') ~ 

( 40 ) 


n{r) - Oqo G+{id')G-{i'd)G+{jG')G-{j'G) - »oo 
r(r) G+{jd)G+{i'G')G-{iG)G-{j'd') ’ 

(41) 

should scale purely logarithmically at the plateau tran¬ 
sition, either as logr or as log^ r. It would be re¬ 
ally interesting to compute this quantity numerically in 
the Chalker-Coddington model and to try to fit it with 
a log^ r + log r + 7 to check this prediction. 

Expanding Green’s functi ons g raphically in terms of 
retarded and advanced path^^^^, these logarithmic cor¬ 
relations appear from the mixing of the operator creat¬ 
ing two advanced paths and two retarded paths (which 
is essentially the analog of a watermelon operator, with 
dimension ^ 2,2 in the notations of Ref. | 66 ]), and the op¬ 
erator which is just counting the number of 

retarded and advanced paths around the neighborhood 
of r^, corresponding to the “thermal”-like perturbation 
driving the system out of criticality by taking z ^1. 


C. Remarks on open vs closed quantum networks 
and conductance correlations 

We emphasize that our field theory predictions ap¬ 
ply to closed quantum networks, for which no lead is 
connected to th e system. In particular, quantities like 
G+(i, i) (diffusion propagator) in terms of the 

Green’s functions ( p! 8 | ) require an infrared regulator in 
order to be finite (an infinitesimal level broadening mak¬ 
ing 1 2 ) I slightly less than 1). Erom a numerical point of 
view, it is more convenient and efficient to work with 
open quantum networks that are also more natural in 
the context of transport. Eor example, to define the point 
contact conductance g{i^j) , one simply cuts in half two 
links i and j (the contacts) with each cut leading to one 
in-going link (source of current) and one out-going link 
(drain). This amounts to attaching two leads to the sys¬ 
tem. The conductance can then be computed within the 
Landauer-Biittiker formalism in terms of the transmis¬ 
sion matrix between the leads (see Ref. for details). 
Unfortunately, opening links in the network correspond 
to the insertion of additional operators that make our 
analysis more complicatecP. Moreover, the point-contact 
conductance g{i,j) is in general completely unrelated to 
G+(i, j)G“(j, i) (except in the SQHE case where their 
average are actually identical up to a constant) - see 
in particular the discussion in Ref. |62l It is nevertheless 
tempting to conjecture that the symmetry considerations 
of this paper extend to g{i,j) as well, even if in general 
it scales in a different way as G+(i, j)G“(j, i). This con¬ 
jecture is motivated by the naive guess that the open or 
closed nature of the network model should not affect our 
symmetry-based analysis: it is then natural to expect 
logarithms as well in observables involving g{i^j) instead 
of G+(i, j)G“(j, i). The analog of eq. ( [38| ) in terms of 
conductances would give an observable characterizing the 
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correlations between the local transport properties (con¬ 
ductances between i and i', and between j and f) of two 
remote regions and rj. It would be very interesting to 
investigate numerically (or theoretically using a different 
approach) whether such physically appealing observables 
show logarithmic correlations at criticality. 


V. DISCUSSION 


In this paper, we have generalized the ideas of Card}!^ 
to analyze the physical origin of logarithmic correlations 
in Quantum Hall plateau transitions using the replica 
trick. In the case of the Spin Quantum Hall transition, we 
found a combination of Green’s functions that scales log¬ 
arithmically on average at the critical point (see eqs. (29) 
and Alternatively, we recovered this result inde¬ 

pendently using the percolation description of the SQHE 
transitiorP^ and a replica-like analysis of the n ^ 1 limit 
of SU{n) dense loop models. Using a similar argument 
for the Integer Quantum Hall transition, we uncovered 
a relatively simple observable that should have logarith¬ 
mic correlations at the plateau transition (see eqs. (38) 
and ([^). 

We emphasize that these predictions based on the 
replica limit of a global symmetry group at the criti¬ 
cal point should be taken with care: as we have already 
mentioned repeatedly throughout this manuscript, con¬ 
sidering a symmetry group smaller than the actual sym¬ 
metry of the critical point may in principle completely 
spoil our predictions: different representations for the 
smaller group that are mixed and lead to logarithms in 
the replica limit n ^ 0 could in fact be part of the same 
larger irreducible representation for the actual symme¬ 
try group of the system. In other words, the amplitudes 
in front of the logarithms within our approach could in 
principle be zero. However, we recall that the very same 
approach lead to results in very good agreement with nu¬ 
merical simulations^ and with more rigorous algebraic 
approaches^ in simpler critical points including the two 
dimensional percolation problem studied in Sec. [Hj even 
though the actual symmetry there is also much larger 


than SU{n) as welP^. Moreover, the replica approach 
seems to lead to very natural predictions for the scal¬ 
ing operators of the critical theory, with results that are 
consistent using two different approaches (replica trick vs 
supersymmetry) in the SQHE case. 


Our predictions should therefore be taken as reason¬ 
able conjectures as to where to look for logarithmic cor¬ 
relations in quantum Hall plateau transitions. It would 
obviously be very important to verify the results of this 
replica analysis numerically, and we expect in particu¬ 
lar that the prediction (41) should be very useful in that 
respect, since it should scale purely logarithmically at 
the critical point, with no power-law dominant contri¬ 
bution. Our field theory predictions apply naturally to 
closed quantum networks, and it would be crucial to see 
to what extent they generalize to open networks that are 
more natural from the point of view of transport, and 
that seem more practical numerically. We note that we 
considered simplified o bservab les for the SQHE and the 
IQHE transitions (see HID 2| above): we leave the de¬ 
tailed analysis of the Sp{2n) or U (2n) symmetry of these 
theories for future work. It would also be very inter¬ 
esting to investigate to what extent the replica and the 
supersymmetry approaches agree for the SQHE transi¬ 
tion: for example, the replica analysis provides a natural 
expression for the “watermelon” operators scaling with 
the /c-hull percolation exponents. Whether these expres¬ 
sions agree with the mapping onto percolation via the 
supersymmetry trick remains unknown. We expect that 
pushing further this symmetry-based replica analysis of 
Quantum Hall transitions should be very helpful to make 
progress towards a deeper understanding of the CETs de¬ 
scribing these critical points. 
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